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The large-scale production of graphene monolayer greatly relies on epitaxial samples which often 
display stress-relaxation features in the form of wrinkles. Wrinkles of graphene on Ir(lll) are 
found to exhibit a fairly well ordered interconnecting network which is characterized by low-energy 
electron microscopy (LEEM). The high degree of quasi-hexagonal network arrangement for the 
graphene aligned to the underlying substrate can be well described as a (non-Poissonian) Voronoi 
partition of a plane. The results obtained strongly suggest that the wrinkle network is frustrated 
at low temperatures, retaining the order inherited from elevated temperatures when the wrinkles 
interconnect in junctions which most often join three wrinkles. Such frustration favors the formation 
of multi-lobed wrinkles which are found in scanning tunneling microscopy (STM) measurements. 

The existence of multiple lobes is explained within a model accounting for the interplay of the 
van der Waals attraction between graphene and iridium and bending energy of the wrinkle. The 
presented study provides new insights into wrinkling of epitaxial graphene and can be exploited to 
further expedite its application. 


I. INTRODUCTION 

Graphene is considered to be a promising material for 
various technological applications due to its supe¬ 

rior physical properties [5]. However, the desired in¬ 
trinsic properties are often downgraded in chemically- 
synthesized graphene, including epitaxial graphene as 
the major large-scale source of the material, due to the 
structural defects in the honeycomb lattice (vacancies, 
pentagon-heptagon defects, grain boundaries) and out- 
of-plane deformations, where the latter can also be in¬ 
duced by the former [6]. Characteristic deformations in 
epitaxial graphene are wrinkles which form during the 
synthesis due to the lateral stress induced by the mis¬ 
match between thermal contraction of graphene and its 
substrate. In the literature, the out-of-plane deforma¬ 
tions of graphene are sometimes also referred to as rip¬ 
ples, ridges or folds, depending on their size and struc¬ 
tural details. In the present work, we do not make this 
distinction and by wrinkles we imply delaminated parts 
of graphene which can be thought of as one-dimensional 
(ID) objects at micron-scales, while forming a wrinkle 
network (WN) at larger scales as we demonstrate in the 
following. Wrinkles are found in graphene on metals 
[MSI, silicon carbide [HHis] as well as various sam¬ 
ples based on transferred graphene [izHini. Besides in 
graphene, wrinkling occurs in other types of layered ma¬ 
terials mi ns and more generally in various systems, 
ranging from semiconductors [23] to soft matter and bi¬ 
ological tissues [24|. 

The prominent role of wrinkles is that they affect the 
transport properties of graphene resulting, for example, 
in tunable transport gaps [25], reduction and anisotropy 
of electrical conductivity nmni and a decrease of ther¬ 


mal conductivity m- Moreover, highly curved surfaces 
of wrinkles exhibit increased chemical reactivity which 
has been predicted earlier m and was confirmed in ex¬ 
periments with oxygen etching [28l [29] and the growth 
of organic semiconductor structures [30]. Wrinkles are 
also important for the process of graphene intercala¬ 
tion. They facilitate the penetration of atoms across the 
graphene sheet either through defective regions charac¬ 
teristic for wrinkles interconnection sites mi ED or via 
temporary defects which may form due to an increased 
reactivity at curved graphene areas [33]. They also serve 
as pipelines for diffusion of already intercalated species 
[29[ |34]. On a macroscopic scale, wrinkles play a sig¬ 
nificant role in the wettability and optical transmittance 
of graphene m- Therefore, a control over their macro¬ 
scopic ordering accompanied by a systematic character¬ 
ization is desirable for a successful technological imple¬ 
mentation of graphene. Some efforts have already been 
made in the direction of wrinkle manipulation, e.g. modi¬ 
fication of wrinkle orientation and length by the scanning 
tunneling microscopy (STM) tip for graphene on silicon 
carbide [15] or alteration of the wrinkle area density and 
orientation by selecting the appropriate substrate mor¬ 
phology for transferred graphene [I2l|36]. 

From the published studies, typical heights and widths 
of wrinkles in epitaxial graphene are found to vary from 
below one and up to several tens of nanometers, de¬ 
pending on the substrate and graphene growth procedure 
[iQi HB da iisi na HE]. Internal substructure of wrin¬ 
kles is often approximated by a single-lobed, semicircular 
shape [iniiisiin]. However, other physical interactions 
involved in the wrinkle formation (e.g. van der Waals) 
may lead to the formation of more complex structures 

[HIST]. 
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Among epitaxially grown systems, graphene on irid¬ 
ium (111) surface [Gr/Ir(lll)] is widely-explored due to 
the high structural quality [38] and almost intrinsic elec¬ 
tronic properties iMiiin] of graphene. Up to now, ex¬ 
perimental characterization of wrinkles of Gr/Ir(lll) via 
low-energy electron microscopy (LEEM) and a simple 
model accounting for the wrinkle formation have been 
reported in the work of N’Diaye et ah |1Q| . There the 
authors focused on graphene islands and graphene con¬ 
taining several rotational domains with respect to the 
substrate. Similar samples were also used in the study of 
Loginova et al. El- Moreover, an x-ray standing wave 
studies revealed correlation between graphene coverage 
and its periodic buckling amplitude mnn] which is a 
consequence of the partial stress relaxation in graphene 
islands via their lateral expansion. Therefore, the struc¬ 
ture of wrinkles is also expected to vary with graphene 
coverage. A detailed study on the relation between lattice 
parameters of graphene and iridium in a broad tempera¬ 
ture interval revealed a hysteretic behavior of wrinkle for¬ 
mation and their flattening [43] . The same study showed 
that wrinkles on average relax ^ 2/3 of the compressive 
stress imposed by iridium while ^ 1/3 still remains in 
the graphene layer. 

In this work, we are interested in wrinkles of the 
high-quality macroscopic graphene samples and hence 
we focus on the full monolayer of Gr/Ir(lll) and refrain 
from the description of graphene islands and incomplete 
graphene layers. LEEM is used for the characterization 
of WN of Gr/Ir(lll). The in situ imaging at elevated 
temperatures gives a detailed insight into the process of 
wrinkling. The WN is found to be well represented by 
an appropriate Voronoi partition of the plane with cen¬ 
troids constructed from the wrinkle positions. Analysis 
of STM measurements reveals that individual wrinkles 
have intriguing and complicated profiles. These are ex¬ 
plained within a framework of a ID model accounting for 
an interplay of different energies important for the wrin¬ 
kle formation. Comparison of our model with experimen¬ 
tal data strongly suggests that the wrinkles we observe 
are multi-lobed, i.e. that the graphene approaches and 
delaminates from iridium several times within the span 
of a single wrinkle. 

II. EXPERIMENTAL 

Sample preparation and all experimental measure¬ 
ments were carried out in situ under ultra-high vacuum 
conditions. The Ir(lll) single-crystal was cleaned by ar¬ 
gon ion sputtering at room temperature at 1.5 keV ion 
energy, followed by oxygen glowing at 850 °C and anneal¬ 
ing at 1200 °C. Graphene synthesis consisted of temper¬ 
ature programmed growth (TPG) followed by a chem¬ 
ical vapor deposition (CVD) which are known to yield 
full coverage of high quality graphene monolayer of uni¬ 


form orientation [33] as confirmed by low-energy electron 
diffraction (LEED) after the synthesis. During TPG, 
ethylene was adsorbed on iridium at room temperature 
(1.33 X 10“^ mbar for 45 s), followed by a flash to 1200 
°C. During CVD, the sample was kept at 850 °C while 
being exposed to ethylene (1.33 x 10“^ mbar for 900 s). 
The LEEM measurements were performed at the U5UA 
beamline of the NSLS synchrotron (Brookhaven National 
Laboratory) with an Elmitec LEEM III system. Sam¬ 
ple temperature (T^), imaging electron energy {Ei) and 
field of view (EOV) varied and are noted in figure cap¬ 
tions. All LEEM images were recorded in the bright field 
mode. Spatial resolution of the microscope was ^ 10 nm. 
STM experiments were performed in Zagreb using Specs 
Aarhus STM operating at room temperature. STM data 
was post-processed in WSxM software |45) . 

III. RESULTS AND DISCUSSION 
A. Wrinkle network 

Eigure [pi shows a LEEM image of an area of 
Gr/Ir(lll) where the dense-packed rows of Ir are aligned 
to the zig-zag directions of graphene lattice, referred to 
as RO graphene [36]. Graphene wrinkles are identified as 
thick dark lines extending ^ 1 /rm in length and forming a 
network. Thinner and more closely spaced wavy lines are 
steps of the iridium surface. It can be noted that wrinkles 
extend mainly in three directions rotated by an angle of 
^ 60°. As a consequence, wrinkle junctions most often 
involve three wrinkles separated by 120° rotation. This is 
directly confirmed in the Eourier transform (ET) of the 
background-subtracted Eigure [^, where a clear 6-fold 
symmetrical pattern is revealed, as shown in Eigure [^. 
The ET pattern consists of stripes and lacks spots which 
indicates that the separation between parallel wrinkles 
and their apparent thickness have a fairly broad distribu¬ 
tion, but the directional preference is nevertheless clear. 
In order to perform a more quantitative analysis, radial 
sums of the ET are presented in a polar graph in Eigure 
[^. Again, the 6-fold symmetry is evident. After plotting 
the same data in a rectangular graph and fitting the six 
peaks with Gaussian curves (not shown), one can pre¬ 
cisely determine the main directions in which the wrin¬ 
kles extend. They are (53.5 ± 0.3) °, (112.5 ±0.3)° and 
(172.6 ±0.6)° (cf. yellow lines in Eigure [^). The dif¬ 
ference between these directions fits well to 60°, which 
agrees with the initial inspection of the LEEM image in 
Eigure [^. This kind of ordering along with the obser¬ 
vation that wrinkle junctions contain mostly three wrin¬ 
kles suggests that the wrinkles form a network composed 
dominantly of hexagons. 

The exact orientation of the dense-packed atomic rows 
of the iridium surface, determined by precise movement 
of the sample along the high-symmetry directions dur- 
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Figure 1. (a) LEEM image of Gr/Ir(lll). Thick dark lines 
extending over the entire sample surface are graphene wrin¬ 
kles, thinner wavy lines are iridium steps. EOV = 9.3 //m, 
Ei — 16 eV, Ts = 96 °C. (b) Eourier transform of image from 
(a). (c) Polar plot of radial sums extracted from (b). The red 
curve is plotted with a negative offset and in arbitrary units, 
(d) Arrangement of iridium atoms (gray spots), graphene (or¬ 
ange hexagons) and dominant directions of graphene wrinkles 
extension (yellow dashed lines). 


ing LEEM imaging, is schematically illustrated in Figure 
[^. In addition, yellow dashed lines mark three deter¬ 
mined wrinkle-extending directions and a perfect coinci¬ 
dence between them and dense-packed rows of iridium is 
found. The shrinkage of iridium during cooldown after 
graphene synthesis can be simply described as a reduc¬ 
tion of separation between adjacent atomic planes in the 
bulk and also atomic rows at the surface. From the obser¬ 
vation of wrinkles directional ordering, it can be inferred 
that the shrinkage induces stress which is relaxed through 
the formation of wrinkles in the direction parallel to the 
atomic rows of iridium i.e. parallel to the zig-zag direc¬ 
tion of graphene. Similar orientational locking was also 
proposed earlier for graphene flakes on iridium m- 

For a comparison to the preferred RO graphene, orien¬ 
tation of the R30 graphene wrinkles was examined where 
the dense-packed rows of iridium are parallel to the arm¬ 
chair directions of graphene. Figure shows the LEEM 
image of a boundary area between the RO and R30 do¬ 
mains of graphene, as identified via micro-LEED patterns 
(upper panels), each exhibiting different LEEM reflectiv¬ 
ity. The underlying Ir steps and wrinkles are present on 
both domains. The corresponding FT patterns from dif¬ 
ferent graphene rotations are shown in bottom panels of 
Figure Whereas the symmetry of the RO transform is 



Figure 2. LEEM image of R0/R30 graphene boundary (cen¬ 
ter) with the corresponding micro-LEED diffraction patterns 
(upper panels) and Eourier transforms (bottom panels). EOV 
= 14 /xm, Ei = 3.8 eV, T, = 98 °C. 

6-fold, the same as in Figure [^, this is not the case for 
the R30 graphene where only two dominant, mutually 
perpendicular directions with a large angular spread are 
found. In principle, two distinct directions are enough for 
stress relaxation in two dimensions and the observation 
of three (two) directions in the RO (R30) graphene calls 
for an explanation. We note that the graphene bend¬ 
ing rigidity is essentially the same along the zig-zag and 
armchair directions in the range of bending radii (up 
to ^ 10 nm) relevant for wrinkling m, ruling out the 
graphene lattice itself as a cause of the observed sym¬ 
metry difference. It was argued in previous studies that 
the R30 graphene is more weakly bound to Ir as com¬ 
pared to the RO graphene [48]. This suggests that (i) the 
RO graphene is compelled to follow contraction of the 
substrate in each of the three equivalent directions de¬ 
fined by the (111) surface of FCC lattice and (ii) the R30 
graphene does not have to keep track of the underlying Ir 
and is free to wrinkle in two approximately perpendicular 
directions. Apparently, the preference between the two 
wrinkling regimes is governed by the binding strengths 
of RO and R30 graphene. However, the absolute orien¬ 
tation of the two wrinkling directions of R30 graphene 
is related to the substrate. It is visible in Figure that 
some of the R30 wrinkles are parallel to Ir steps while 
others are perpendicular to them. Similar observations 
were made in the case of graphene on Pt(lll), a system 
in which graphene is more weakly bound to the substrate 
as compared to the RO Gr/Ir(lll) [71 [29], 


B. Voronoi diagram approximation 

On the basis of collected LEEM data, we can show 
that the WN on the RO graphene can to a good extent 
be approximated by a Voronoi diagram (VD) i.e. a net¬ 
work obtained in a Voronoi tessellation of a plane |49| . 
The construction of a VD starts from a set of points 
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Figure 3. LEEM image of the graphene WN superimposed 
by the VD (yellow dashed lines) generated by red centroids. 
Arrows indicate several deviations of the diagram from the 
network. The inset shows Eourier transform of the centroids. 
Image size: (6.5 x 4.3) /xm^, Ei = 1.9 eV, Tg = 91 °C. 


(called centroids) which divide the plane into cells (called 
Voronoi cells) associated to every centroid in a way that 
each cell contains only points that are closer to a given 
centroid than to any other. All Voronoi cells together 
constitute a VD. Figureshows the recorded LEEM im¬ 
age of Gr/Ir(lll) overlapped with the VD (yellow dashed 
lines) generated from a set of centroids marked as red 
dots. The positions of centroids can be determined from 
the WN by a simple geometrical construction applied to 
each cell of the network m (see Supplementary data for 
details). In the case of WN, large number of wrinkles 
coincide with the Voronoi cell borders with some minor 
exceptions. For example, the VD predicts the existence 
of wrinkles which are actually not there (cf. arrows 1 and 
2) or the absence of wrinkles which are found at certain 
locations (cf. arrows 3 and 4). In spite of this, the VD 
approximation is satisfactory and is likely to get even bet¬ 
ter in cases when various defects on the iridium surface 
as well as iridium steps are suppressed to a minimum. 
In addition, shorter and smaller wrinkles are hard to re¬ 
solve in LEEM and wrinkles in general can be mistaken 
for iridium steps and vice versa, which may deteriorate 
the approximation to some small extent. 

It might seem at first that the spatial distribution of 
centroids in Figure is random, but if that were the 
case wrinkles would not exhibit preference in their direc¬ 
tion and the corresponding Fourier transforms would be 
isotropic, which is in contrast to what we observe. More¬ 
over, Fourier transform of the centroids alone, displayed 
in the inset of Figure exhibits hexagonal symmetry 
and indicates their quasi-periodic arrangement. In fact, 
a typical Voronoi cell of a WN is hexagonal, as can also be 
concluded from the analysis of Figure It follows that 
the WN as well as the corresponding VD are networks 
of distorted hexagons similar to a honeycomb structure. 


Based on these findings, WN can be simulated via con¬ 
struction of perfectly arranged honeycomb VD to which 
a specific amount of disorder has been introduced (see 
Figure S2 in the Supplementary data). VDs generated 
from a random distribution of centroids are referred to 
as Poisson-Voronoi diagrams and have been studied ex¬ 
tensively IlSEI]- Expected values of various parameters 
(cell area, cell perimeter, number of edges etc.) of such 
diagrams have been calculated and simulated. These val¬ 
ues are not expected to vary significantly for centroid dis¬ 
tributions other than random m, so they can be adopted 
for the characterization of the WN found for wrinkles in 
the RO graphene. It is worth noting that the expected 
value of an angle at a typical wrinkle junction is 120°, 
as we also inferred from the data presented in Figure 
One of the easiest quantities to obtain from LEEM is 
the number of wrinkle junctions per unit area, rij. Once 
this number is known, expected values of several param¬ 
eters of the network can be calculated from the theory of 
Poisson-Voronoi tessellation [49]. For the LEEM image 
in Figure]^ one finds rij = 1.4 /xm'^ and other derived 
parameters are given in Table [I] For comparison, we also 
list the values of other parameters obtained from the ex¬ 
periment (Figure]^ which show good agreement with the 
theoretical predictions. Hence, by identifying the WN as 
a VD, one can easily estimate various parameters rele¬ 
vant for network’s general description which in turn fa¬ 
cilitates evaluation of the effect the wrinkles may have 
on the properties of graphene, e.g. electronic transport. 
Finally, the values given in the table are not universal 
but are expected to differ depending on the preparation 
parameters of graphene, especially the growth tempera¬ 
ture. 


C. Dynamics of wrinkling 

After the formation of 1 ML graphene and cooldown to 
room temperature, the sample was re-heated to 650 °C 
and then slowly cooled down while LEEM images were 
being simultaneously recorded. Such experiment enabled 
the observation of disappearance and reappearance of in¬ 
dividual wrinkles. In Figure]^ and b we show two areas 
of the sample recorded right before (frame 1) and af¬ 
ter (frame 2) the appearance of several wrinkles during 
cooldown. By comparing frames 1 and 2, we deduce that 
the reflectivity of graphene surrounding the appearing 
wrinkles increases by 3% once the wrinkle is formed (see 
also the movie clip in the Supplementary data). This is 
even clearer in frame 3 where the difference of frames 2 
and 1 is shown. This effect is discussed in the work of 
N’Diaye et al. m to be a result of a change in interfer¬ 
ence conditions of incoming and reflected electron waves. 
However, in the following we propose an additional mech¬ 
anism for this local change of reflectivity. 

In our experiment, LEEM images are recorded in 
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Table L Parametrization of the WN. Parameters of the Poisson-Voronoi diagram (1®^ column), their mutual dependence (2^^^ 
column), values extracted from the LEEM image shown in Eigure column) and values calculated from the equations 

given in the second column (4^^ column). Calculations were performed based on the experimentally determined value rij = 1.4 


parameter 

relation to nj 

experimental value 

calculated value 

density of junctions, nj 

rij 

1.4 /im“^ 

- 

density of centroids 

Uj 

2 

- 

0.7 

wrinkle length per unit area 

9 /T7 

2 

1.4 

1.7 Atm“^ 

area of a cell 

2 

rij 

- 

1.5 iiw? 

perimeter of a cell 

V '^3 

- 

4.8 /xm 

length of a single wrinkle 

2 AT 

3 Y ^3 

(0.7 ± 0.3) /xm 

0.8 /xm 



Eigure 4. (a, b) LEEM micrographs, (3x3) fivo? in size, show¬ 
ing two different areas in of Gr/Ir(lll) right before (frame 1) 
and after (frame 2) wrinkle emergence. Erame 3 displays the 
difference 2 — 1. (a) is recorded at 350 °C, (b) at 270 °C. 
Ei — 3.8 eV. (c) An illustration of the origin of LEEM reflec¬ 
tivity changes upon wrinkle formation (see text for details). 
White spots represent (0,0) and moire diffraction spots. Yel¬ 
low circle indicates size of the contrast aperture used for 
LEEM imaging, km is the moire wave vector, its variation at 
temperatures Ti and T 2 is exaggerated for clarity. 


bright field mode at 3.8 eV electron energy and at such 
low energies only the (0,0) spot and a fraction of higher 
order moire spots contribute to image formation (cf. Fig¬ 
ure]^). Furthermore, at the onset of wrinkle formation, 
the surrounding graphene lattice relaxes, leading to an 
increase of the graphene lattice parameter and conse¬ 
quently a decrease of the moire wave vector module km- 
This brings moire diffraction spots closer to the (0,0) spot 
by roughly 3% [43] (cf. Figure]^). Hence, a decrease of 
km results in an increase of the fraction of moire spots in¬ 


tensity contributing to LEEM images. In other words, ar¬ 
eas in LEEM images which become brighter after a wrin¬ 
kle forms are areas where the graphene lattice is relaxed 
as compared to its surroundings. In this way graphene 
lattice relaxation can be directly observed in real space. 
This also indicates that the stress imposed on graphene 
is visibly relaxed only in the ^ 0.5 /rm vicinity of the 
wrinkle and that graphene areas situated further away 
remain highly stressed, at least at elevated temperatures 
of several hundred degrees Celsius. The scenario of local¬ 
ized stress relaxation is supported by recent theoretical 
simulations showing that graphene stretching energy is 
reduced in the wrinkles’ surroundings [52]. Upon further 
sample cooldown, the reflectivity (and hence graphene 
stress) in the entire LEEM field of view equalizes but 
still exhibits local variations which appear to be related 
to iridium surface steps. Similar relation between stress 
non-uniformity and substrate morphology has been ob¬ 
served for graphene on silicon carbide m- Finally, at 
room temperature when the surface is characterized by 
a fully formed WN, the reflectivity variations are dimin¬ 
ished and indicate uniformity of stress in graphene. 

Monitoring of the disappearance and reappearance of 
wrinkles in LEEM during re-heating enabled examination 
of wrinkle nucleation sites. Figures [^a) and (b) show two 
LEEM images of the same area on the sample. Image (a) 
was recorded at 247 °C during heating and image (b) at 
247 °C during cooling of the sample. In between, the 
sample temperature reached a maximum of 650 °C when 
all wrinkles were flattened out. Re-heating cycles involve 
thermal drift of the sample, but prominent iridium fea¬ 
tures (such as a series of distinct circular terraces, cf. red 
arrow in Figure enabled identification and observation 
of the same surface area for the given temperature. Com¬ 
parison of the two images allows us to infer that some of 
the wrinkles reappeared at the same location (cf. white 
arrows in Figure but many other formed at new posi¬ 
tions which were wrinkle-free before re-heating (some of 
them marked by yellow arrows in Figure]^. The num¬ 
ber of wrinkles in Figure is apparently lower than in 
Figure which is due to the hysteretic behavior of the 
wrinkle formation [43]. Additional wrinkles appear as the 
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Figure 5. (a, b) LEEM images illustrating wrinkle pinning. Several graphene wrinkles are visible (a) at 247 °C during heating 
and (b) at 247 °C during cooling. Some of the wrinkles reappeared in the same locations (white arrows) while some formed in 
new locations (yellow arrows). Characteristic Ir terraces (red arrows) facilitate identification of the same surface area of the 
sample. EOV = 7.2 /xm, Ei — 3.8 eV. (c) STM image of the wrinkle junction anchored at leftover dirt particle (first derivative 
of the topography image). Tunneling voltage: 630 mV, tunneling current: 1.7 nA. 


sample is further cooled down to room temperature, but 
due to mentioned thermal drift a more precise compari¬ 
son for exactly the same surface area cannot be made. We 
also observe that new winkles often form at the endpoints 
of already existing ones and the whole process sometimes 
resembles to ’’crack” (wrinkle) propagation (see also the 
movie clip in the Supplementary data). 

We conclude from our experiments that various struc¬ 
tural features of graphene or iridium can anchor some 
of the wrinkles, serving as nucleation sites for their for¬ 
mation. However, many wrinkles are exempt form this. 
A typical wrinkle junction containing dirt particle, most 
likely a leftover from the crystal cleaning process, was 
characterized with STM and is shown in Figure [^. Such 
sites are suitable for the nucleation of wrinkles due to 
the reduced binding of graphene to the substrate and 
a presumably large number of defects contained in the 
graphene lattice itself. If dirt particles are stable at el¬ 
evated temperatures, they can anchor wrinkles i.e. they 
are preferred sites of wrinkle nucleation. This is in line 
with previous findings on graphene flakes where the au¬ 
thors claim that majority of wrinkles nucleate at defects, 
more precisely at heptagon-pentagon pairs m- 


D. Multi-lobed wrinkles 

We now proceed to examine the features of individ¬ 
ual wrinkles, in particular their widths and heights and 
the variation of graphene-iridium separation within a sin¬ 
gle wrinkle, i.e. the wrinkle profile. A simple, single- 
lobed profile is often considered, but one may exam¬ 
ine more complicated profiles consisting of several con¬ 
secutive lobes. Motivation for such profiles is found in 
the STM measurements (shown in Figure]^ and below, 
in Figure and it is of interest to study the energet¬ 
ics of their formation. The energies of different pro¬ 
files are compared using a simple ansatz for the pro¬ 
file shape which has been shown to accurately describe 
the compression-governed delamination of graphene [54] 


and other thin films [55l l56] . The profile is represented 
by the cosine function of the amplitude Am extending 
over m periods of length Im^ leading thus to m successive 
identical lobes. The total width of the profile is li (the 
width of a single-lobed wrinkle) for every m. The out- 
of-plane displacement of m-lobed wrinkle in the interval 
X G [—/i/2,/i/2] can be expressed as 


(a;) = ^ 


^ , 27rx 

1 -h cos I — - qn 


( 1 ) 


where g' is 0 if m is odd and is 1 if m is even. An example 
is given in Figure]^ where besides a basic, single-lobed 
wrinkle (m = 1, black line) we also show wrinkles hav¬ 
ing two (m = 2, red line) and three (m = 3, blue line) 
lobes. In order to be able to directly compare wrinkles 
containing different numbers of lobes, the total length of 
wrinkled graphene in the model is kept constant. This 
condition in turn yields the periods and amplitudes of 
multi-lobed wrinkles: Im = h/m and Am = Aijm (cf. 
Figure 1^, see Supplementary data for details). 


As the total delaminated graphene length in the model 
is the same in all profiles, the strain within graphene and 
the accompanying strain energy do not change for differ¬ 
ent profiles. Thus, the profitability of multiple lobes de¬ 
pends only on the balance between the unfavorable bend¬ 
ing energy cost and the favorable van der Waals binding 
energy gained in the areas where graphene is brought 
close to iridium. In principle, additional energy can also 
be gained from the van der Waals interaction between 
the neighboring lobes m- However, for typical profiles 
considered here, this contribution in negligible and is not 
taken into account. We model the van der Waals energy 
of the profile by summing the Lenard-Jones 6-12 inter¬ 
actions between the atoms in the graphene and iridium, 
so that the binding energy per unit area of the profile is 
given by 
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(a) 




Figure 6. (a) Cross-sections of multi-lobed wrinkles for m = 1,2 and 3 (see text for details). Hq is the equilibrium graphene- 
substrate separation, (b) The sum of binding and bending energies of a wrinkle per unit area as a function of number of 
lobes within a wrinkle. Energy minima are found for different number of lobes, as indicated by black dots. The numbers in 
brackets correspond to = (rnAm^Tnlm) in nanometers and optimal m is determined for each curve from the position 

of its minimum. The red dashed curve is plotted in accordance with the parameters obtained from the STM data shown in 
Figure 


EvdW = — 



hp 

Wm{x) + hp 


3 


1 / hp 

2 \Wm{x) + hp 



( 2 ) 


where hp is the equilibrium graphene-iridium separation 
which is set to 3.4 A mi, and 7 is the graphene-iridium 
(111) binding energy per unit area of a flat graphene; it 
amounts to 50 meV per C atom (0.308 J/m^) [41]. The 
bending energy per unit area is given by [58] 



-b/2 


where D is the graphene bending rigidity having a value 
of 0.238 uN-nm [471 EH ETj- The comparison of the to¬ 
tal energy Ey^tw + ^6 for different number of lobes m is 
shown in Figure for several combination of the param¬ 
eters (Mi,/i) = (additional details of the 

calculation are given in the Supplementary data). Black 
dots mark the minima of the displayed curves. The graph 
shows that for certain combinations of li and Ai it is en¬ 
ergetically more favorable to form multi-lobed wrinkles. 
For example, a single-lobed wrinkle with Mi = 5 nm 
and /i = 20 nm [cf. (5, 20) curve in Figure [^] is less 
preferable than a wrinkle containing the same amount of 
graphene but having two lobes so that its amplitude and 
wavelength are A 2 = Ai/2 = 2.5 nm and /2 = ^i/2 = 10 
nm, respectively. In general, multi-lobed wrinkles be¬ 
come preferable when the width of a wrinkle is much 
larger than its amplitude. In such cases, energy gained 
by additional binding of graphene to iridium upon forma¬ 


tion of several lobes compensates for the cost of graphene 
bending. 

The multi-lobed wrinkle structures are observed for 
Gr/Ir(lll). Figure IT shows a constant-current topo¬ 
graphic STM image of a wrinkle, indicating the existence 
of several lobes. In order to resolve the structure in more 
detail, first derivative of the topographic image is shown 
in Figur^^. Additionally, Figure shows the wrinkle 
profile (taken along the red line in Figure]^) where we 
identify four separate lobes as schematically indicated 
by a dashed line. In our STM experiments such complex 
wrinkle profile was often found. Possible multiple-tip ef¬ 
fects are eliminated on the basis of the observation that 
multi-lobed structures are observed in joining wrinkles 
protruding in different directions (cf. Figure Iff) , thus 
regardless of the scanning direction. Moreover, images 
of surrounding fiat areas do not give any indication of 
multiple-tip imaging. The measured topographic wrin¬ 
kle profile does not match our model from Figure [^ide¬ 
ally. One should take into account, however, that our 
model profile is a variational ansatz used only to demon¬ 
strate the propensity towards formation of multi-lobed 
wrinkle profiles. As such, it is not an exact solution of 
the problem, and its comparison with the experiments 
should be performed with care. In addition, the finite 
radius of the STM tip may not allow for clear imaging 
of narrow trenches in between the individual lobes. Fur¬ 
thermore, during STM imaging of such curved structures 


































d (nm) 

Figure 7. STM measurements of multi-lobed wrinkle. Topog¬ 
raphy is shown in (a) and the corresponding first derivative 
is shown in (b). Panel (c) shows a profile across the red line 
in (a) where four lobes are identified as schematically indi¬ 
cated by black dotted line. Blue arrows mark total width and 
average height of the wrinkle. Tunneling voltage: 630 mV, 
tunneling current: 2 nA. 

an increased contributions to the tunneling current from 
the side of the STM tip can be expected, which can af¬ 
fect the measured topographic profile. One should note, 
however, that the experimentally determined lobes have 
similar widths, in reasonable agreement with the model 
calculation which assumes the same widths for all the 
lobes. 

Total width of the wrinkle is 14.36 nm whereas the 
height is estimated to be 0.59 nm on the basis of its 
average out-of-plane displacement (cf. blue arrows in 
Figure [^). In Figure [^, the red dashed line corre¬ 
sponds to these experimental values: li = 14.36 nm and 
Ai = 4 • A 4 = 2.35 nm. Minimum of the curve is found 
at m 3, which is in rough agreement with four lobes 
inferred from the STM data. At this point we would 
like to highlight two important virtues of our model: (i) 
it predicts multi-lobed wrinkles, and (ii) it does so in a 
semi-quantitative fashion, as can be seen from the antic¬ 
ipated number of lobes for the wrinkle shown in Figure 

m 

An issue which remains to be explained is the ap¬ 
pearance of several lobes in a profile of a single wrinkle. 
Namely, lobes could be conceptually separated from each 
other forming thus several single-lobed wrinkles which 
would have the same energy as a multi-lobed wrinkle. 
However, whereas this would indeed be the case for par¬ 
allel wrinkles in essentially ID model considered here. 


the inclusion of separated lobes in the already formed, 
connected two-dimensional WN would require additional 
energy to form new wrinkle junctions. The fact that we 
do observe multi-lobed wrinkles suggests also the sce¬ 
nario for their formation. As the temperature is lowered, 
at a transition temperature the strain in the graphene re¬ 
laxes through the formation of, likely, single-lobed wrin¬ 
kles. The transition from a flat graphene to a WN is not 
instantaneous (see the movie clip in the Supplementary 
data), as there is disorder in the system (it is thus dif¬ 
ficult to speak of a well-defined critical temperature for 
wrinkle formation). As the temperature is further low¬ 
ered, the existing wrinkles grow, the new ones form, and 
they soon connect to each other forming the WN with 
many junctions. Upon further temperature decrease, ad¬ 
ditional strain accumulates. Its relaxation may proceed 
either by formation of new wrinkles, by repositioning of 
the existing wrinkles or by influx of graphene in the al¬ 
ready formed WN, increasing thus both wrinkle ampli¬ 
tude and width. Although a different network topology 
may be energetically preferred at lower temperatures, its 
formation would require a global reconstruction of the 
established WN, so that the additional strain relaxes in 
the already formed network. The dynamical scenario 
that we propose indicates that the topological proper¬ 
ties of the network are fixed in some interval of elevated 
temperature required to form most of the wrinkle junc¬ 
tions. The once formed network frustrates the system 
at lower temperatures, the wrinkle junctions behaving as 
rigid anchors, and the strain relaxes only by increase of 
the width and amplitude of the wrinkle profiles in the 
network established at elevated temperature. The influx 
of graphene in the wrinkle profiles with the lowering of 
the temperature induces the formation of multiple lobes. 
The frustration of the WN is also consistent with the 
observed hysteretic behavior in the heating and cooling 
cycles. 

Even though they were not discussed in the existing lit¬ 
erature in detail, wrinkles exhibiting multiple lobes were 
also observed in transferred graphene grown on Cu foils 
m- In that case, wrinkles are found in a broad range 
of widths. For wide wrinkles, wrinkle height is approx¬ 
imately constant and relatively low (^ I nm), whereas 
for narrow wrinkles it is several times larger (^2 — 6 
nm). This observation can be explained with the aid of 
our model: wide wrinkles tend to form multi-lobed struc¬ 
tures (wider wrinkle ^ larger number of lobes) which is 
reflected in their constant height. Moreover, the accom¬ 
panying AFM measurements of wrinkle profiles provide 
strong indication of a presence of several sequential lobes 
within wide wrinkles. The concept of competition be¬ 
tween binding and bending energies in wrinkling phe¬ 
nomena is also valid for other layered van der Waals 
materials. It was recently shown that the out-of-plane 
deviations found in M 0 S 2 layers exhibit configurations 
similar to multi-lobed wrinkles which are stable only if 
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the total width of such structures is sufficiently large [H] . 
In the case when separation between individual lobes is 
reduced, they merge and form a single-lobed wrinkle due 
to the lower cost of bending energy. 


IV. CONCLUSION 

For the RO graphene, the WN exhibits quasi-hexagonal 
ordering with wrinkles extending in three preferential 
directions aligned with the dense-packed atomic rows 
of the Ir(lll) surface. Such ordering is absent in the 
R30 graphene due to the reduced graphene-iridium bind¬ 
ing. We find that the WN can be well approximated 
by the VD which greatly facilitates its description and 
parametrization. The formation of wrinkles at elevated 
temperatures is accompanied by an increase of LEEM re¬ 
flectivity in their immediate vicinity and we propose that 
the reason for this is localized, inhomogeneous stress re¬ 
laxation of the graphene lattice. Moreover, structural 
features of iridium or graphene, such as wrinkles them¬ 
selves or dirt particles, can act as nucleation sites for 
other wrinkles. Einally, our STM results and model cal¬ 
culation indicate that multi-lobed wrinkle profiles are en¬ 
ergetically stable and result from the frustration of the 
WN during sample cooldown. Our study brings new in¬ 
sights into the properties of wrinkles and WN in epitaxial 
graphene and is of potential use in special electronic, opti¬ 
cal or mechanical applications of graphene and follow-up 
materials. 
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APPENDIX A. SUPPLEMENTARY DATA 

Supplementary data associated with this ar¬ 
ticle can be found, in the online version, at 
http://dx.doi.org/10.1016/j.carbon.2015.07.059. 
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